how to solve 2d heat conduction